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The effective potential of scalar QED is computed analytically up to two loops in the Landau 
gauge. The result is given in 4 — e dimensions using minimal subtraction and e-expansions. In three 
dimensions (e = 1), our calculation is intended to help throw light on unsolved problems of the 
superconducting phase transition, where critical exponents and the position of the tricritical point 
f"") ' have not yet found a satisfactory explanation within the renormalization group approach. 
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$_l 1 I. INTRODUCTION 

a.: 

Recently, important progress has been made to the evaluation of Feynman diagrams in 4 — e dimensions. In 
\Q • particular, two-loop Feynman diagrams with unequal masses of the internal lines are now available analytically 

Moreover, the full e-expansion of the sunset diagram is known This is important for critical phenomena in 4 — e 
dimensions since, for instance, a three-loop calculation requires the knowledge of all terms of order e of the two-loop 
diagrams. 

With the help of these results, we investigate in the present paper an 0(7V)-symmetric version of scalar QED 
containing N/2 complex fields in 4 — e dimensions, which, for N = 2, reduces to the Ginzburg-Landau model of 
superconductivity. Our results are intended to shed some light on the open problems of the superconducting phase 
■ transition in three dimensions. From a dual formulation of the Ginzburg-Landau theory we know that there must 
exist a regime of the parameter k, the ratio between scalar and vector masses, where the transition changes from first 
to second order H. This scenario has been confirmed by several Monte-Carlo simulations, see 

In contrast to this, the traditional 4 — e -dimensional renormalization group analysis of the Ginzburg-Landau 
model has produced more puzzles than answers. At the one-loop order, a second-order transition is obtained only 
for N > N c = 366, much larger than the N = 2 -value of a superconductor HQ. The situation does not improve by 
going to two loops Q . Only a somewhat contrived resummation procedure leads to a desired fixed point |J . Only in 
direct three-dimensional calculations has a second-order transition been obtained for N — 2, see for example |lCj. The 
absence of a charged fixed point for N — 2 in 4 — s dimensions seems thus to be a specific weakness of this approach. 

At present it is hoped that variational perturbation theory may be able to locate a fixed point of the Ginzburg- 
Landau model, thus allowing to extract physical values independently of Ref. JG]. This theory developed in Rcf. Jfi 1 ] 
has proven to be a powerful tool for determining critical exponents in three [^2|-fi4 as well as in 4— e dimension s |1l5| , |l6f 
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Recently, we have applied it to the determination of amplitude ratios in three dimensions of the 0(iV)-model 17 . We 
did these calculations using a method proposed by the Aachen group fl8[| , where analytic renormalization is applied in 
the form of minimal subtraction although working directly in D — 3 dimensions, without any e-expansion. A similar 
determination of amplitude ratios in 4 — e dimensions with e-expansion is still missing. For determining the associated 
amplitude ratios, we need the effective potential in 4 — e dimensions. This is what we want to calculate here. 

As a cross check we shall often set the electric charge equal to zero and check that we recover the results for the 
usual 0(-/V)-symmetric </> 4 -theory. 

The renormalization constants and various beta functions of scalar QED have been calculated up to two loops in 
the early work on this subject [p|,|l9||. It is not our purpose to rederive these results. However, in the process of 
obtaining the effective potential, i.e., of calculating the various Feynman diagrams in the symmetry-broken phase, we 
obtain not only the effective potential but also the renormalization constants of mass Z m i and (/> 4 -coupling constant 
Z g to up to two loops. The renormalization constants of the scalar and vector fields Z$ and Za, on the other hand, 
enters only with the first loop order. The fact that the charge coupling renormalization constant Z e is not obtained 
to the same order as the <fi -coupling constant one has its origin in the fact that, due to a Ward identity, Z e = Zq, 
see Eq. @. 

Before closing the introduction, we point out that the effective potential of the electroweak part of the standard 
model has already been calculated to the two-loop order |IJ. However, our work is more than the abelian U(l) subset 
of this reference for two reasons: first, we work with A72-complex scalar fields (N=2 in |l]]); second, in view of the 
application to critical phenomena in three dimensions, the e-expansion is used. At the two-loop order, this means 
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that the one-loop part of the effective potential is expanded up to the order s. This term is of course not there when 
being interested in going in four dimensions, as is the case ofRef. §. 

Finally, we mention that our work can be seen as the extension to a two coupling constants problem of the work 
of Brezin et al. |2(J (see also which consider the e-expansion of the equation of state of the iV-components (p 4 

theory to the two-loop ordeif]. Without the gauge field, the effective potential can be used to determine this equation 
of state. 

The paper is organised as follow: in Section^, we specify the model and our conventions. In Section [II, we present 



the various intermediate steps for the obtaining of the effective potential. The individual results are combined in 
Section IV, and the conclusions drawn in Section (v[ 



II. MODEL 



The Lagrangian density to be studied contains N/2 complex scalar fields </> coupled to the abelian fields and 
reads, with a covariant gauge fixing, 

C = \D^\ 2 + mV + ||0| 4 + \F% + ^(dM 2 , (1) 

where = — ieA^ is the covariant derivative, is the usual field-strength tensor and a is a gauge parameter. 
There is no need to introduce ghost fields to guarantee gauge invariance since these decouple from the system in an 
abelian gauge theory. 

The effective potential will be obtained using the so-called background- field method of DeWitt |2j| . We shift the 
scalar field by an unknown constant <&: <j) ~~ * ^ + 4 1 - This generates new vertices. To simplify the calculation, we shall 
use throughout the Landau gauge a — > 0. This reduces the number of Feynman diagrams and, since a = enforces 
d^A^ = at the Lagrangian level, removes a possible mixing of A^cf) and A^ft terms, thus decoupling scalar and 
gauge propagators. It is further advantageous to use real fields, defining 

0=J_(0! $ = _L($ 1+Z $ 2 ). (2) 
Then the Lagrangian has the expansion around the background field: 
C = Co + \^[G T P T + G L P L \<$> + ±A„D T P^A U 

+ |[(0 2 ) 2 + 4<&<M0 2 )] + e 2 A 2 $0 + \e 2 Al<f + eA^fod^ - <^<fe), (3) 

where <p and $ are now N components real fields written as two-dimensional iso- vectors <p = (</>i, fo) and $ = $2)- 
The notation is: 

(4) 

a 2 + m 2 + | ( j ) 2 5 G L = _Q2 + m 2 L = _ d 2 +m 2 + g$2 j (5) 

d 2 + e 2 * 2 , (6) 



Co 


1 2 
= -m 


$ 2 + i-$ 4 




4! 


G T 


= -d 2 


+ TOy = — 


D T 


= -d 2 


+ m 2 = -< 


P T 

y 


= s a - 


$2 ' 



pl = Mi pt 6 Mil (7) 



with G T , G L being the transverse and longitudinal inverse propagators of the scalar field, and the inverse trans- 
verse propagator of the photon field. The transversality of the latter is due to the Landau gauge. Note that there is 
no term eA M d M (<I>20i — &i4>2) which would mix vector and scalar propagators. 

Compared to the complex-field representation, real fields have one small complication: the gauge-scalar-scalar vertex 
involves a vector product since it mixes the real and imaginary parts </>i and <f>2- For complex fields it is diagonal in 
the scalar field space. 



1 For the Ising model, corresponding to TV = 1, the authors of have succeeded going to the three-loop order of the 
e-expansion of the equation of state. 
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Up to two loops, the fastest way of determining the Feynman diagrams with their proper weight is to use the Wick 
theorem. In higher loops, it might be interesting to use a more efficient algorithm to generate the diagrams and their 
weight, in the same spirit as, for instance, pi 25 1. Using a similar procedure, we have derived the diagrams of scalar 
QED up to the fourth order. Their calculation will form the subject of a separate publication p6fl . 

We shall calculate the diagrams in the following expansion, in which h counts the number of loops: 



^eff. pot. 



^classical 



(w - 1) + +(D - 1) 0] 



(N 2 



uoo^-^oo^ool^-^-^oo+oo] 



(iv-i) <3> + 3 <3 



4 I - \ <I» J 



<t^j > ~ e2 ( N ~ 2 ) + 2 J j. 



(8) 



Explicitly, this equation reads in a unit D-dimensional volume: 



V = -m 2 $ 2 
2 



£$ 4 + -[(N- l)Trln(G T ) + Trln(G L ) + (D - l)Trln(C T )] 



\f — 1) J dp A(p, mr) + J dp A(p, oil) 
J dp A{p, m 7 ) 



+ 2(JV-1) 



dp A(p, ttit) 



dpA(p, m L ) 



(N-l) J d P A(p,m T ) + J dpA(p,m L ) 
-(^) 2 $ 2 (N-l) J dpdqA(p,m T )A(q,m T )A(p + q,m L ) + 3 / dpdqA(p,m L )A(q,m L )A{p + q,m L ) 

dpdq P»v(p) p Jv(<1) A (Pi m J )A(q,m- / )A(p+ q, m L ) 
J dpdq P* v (p) qiM q v A(p, m 7 ) [(N ~ 2)A(q, m T )A(p + q, m T ) + 2A(q, m L )A(p + q, m T )\ J j , 



(9) 



where dp = d D p/(2ir) D , and A(p,m) — l/(p 2 + m 2 ). 

All quantities in these expressions (fields, coupling constants, and masses) are bare quantities. Up to the second 
order in the loop expansion, the divergences show up as poles in e up to the order 1/e 2 . They have to be removed 
to have a finite limit e — > 0. They may either be removed by adding counter-terms to the initial Lagrangian, and 
calculating the associated Feynman diagrams in addition to those in (||). Alternatively, and this is how we shall 
proceed, we may use (@) and (^) and include renormalization constants which are reexpanded up to any given order 
in the loop parameter, i.e., up to the order h 2 in our case. These renormalization constants are defined as 



Zj>, Afj, = A /lr V Z A , 771 = 777 



, Z„ 



z& 



e 9 ell e e >* 

.9 = 5r/J^2, e = e r fi 1 = —= 



e/2 



(10) 



where, in the last equation, we have taken into account the relation Z e — Z$, which is a consequence of a Ward identity. 
Intuitively, it comes from the requirement D^cf) ^fZ^D^r for the covariant derivative. In the above equations, 
the bare quantities are on the left-hand-side and the renormalized ones are on the right-hand-side, indicated by the 
subscript "r" . 

We also note that the vacuum requires a special treatment J27j]. For <I> = 0, the dimensionality requires V cx m D . 
We therefore add a term m A h/g to the Lagrangian. With this we define a new renormalization constant absorbing 
the vacuum divergencies by 

4 4 

(11) 
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In the next Section, we calculate the different loop orders, based on 
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III. EVALUATION OF THE DIAGRAMS 

In the following, the subscript indicating the renormalizcd quantities will be omitted, for brevity of notation. The 
different renormalization constants are expanded with respect to the loop-parameter: 



3 ' 



(12) 



z=i 



The results are then identified order by oder in h. To fix the ideas, we give here to procedure for the gauge field 
trace-In: The bare diagram has to be modified according to: e 2 $ 2 — > e 2 § 2 Z^/ Za, and the result reads, to the first 
order in h 



Trln(D T ) = Trln ( -d 2 + e 2 $ 2 



Z A 



Trln(-3 2 + e 2 $ 2 ) + he 2 &(Z^ ] - Z^)Tr 



1 



-d 2 + e 2 $ 2 



The last term will contribute to the two-loop result, removing parts of its e-poles. 



(13) 



A. Renormalized zero- and one-loop order 

The renormalized zero-loop order is trivial 



^ = 0) = -?7l 2 $ 2 + fr//$ 4 



m 

4!^ + ^' 

while the renormalized one is simply a combination of the previous bare zero-order and the trace-In terms: 



2 ™ z 41' 

r(l - D/2) 1 
(4tt) b / 2 D 1 



,.v , )( <) 2 (4) e/2 +(™i) 2 f4) £/2 +^-i)K) 2 f4 v 



(14) 



(15) 



The constants Z ■ are determined in order to remove the e-poles at D = 4 in the Euler T(l — D/2) function. In four 



dimensions, it would be sufficient to make the e-expansion of (15) up to the order e . However, we are also interested 
in going to dimension D — 3. For this reason, it will be necessary to perform the e-expansion up to the order e. 



B. Renormalized two- loop order 



The two-loop order contains contributions from the zero- and one-loop bare diagrams, which will cancel the two-loop 
poles from the bare two-loop diagrams. These contributions read ("ct" is for counter-terms) 

V(l = 2) ct = -m 2 ^> 2 Z {1 l + ^if&Z^ + —ZW 



r(i-p/2) i _ £ 

(4n) D / 2 2' 
( u 2 \ 6/2 



»-4(N l)(m|) ^ E/2 [±(3*4 - m|)Z« + \{m 2 L - m 2 T )Z^ - m%zf 



-,/) - W) [^^izf-zf) 



(16) 



To simplify the evaluation of the two-loop diagrams, we introduce the functions J(m) and I (mi, rri2, m^), where 
we have kept the same notation as in |5J]: 
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J(m) — j dpA(p,m), 
/(TOi,m 2 ,m 3 ) = / dpdq A(p,mi)A(q,m 2 )A(p + q,m 3 ). 



(17) 
(18) 

The function J(m) is trivial to determine, and its result has in fact been used in jig ). For I (mi, 7712,777,3), we use 
the result obtained in which are in accordance with |g] and the recent work |3j] . The latter reference gives an 
expansion for all order in e but, since we need in the present work only the term of order e , we stop at this order: 



J(mi 

(^ e ) 2 (47r) 4 /(mi,m 2 ,m3 
1 



r(i-i>/2) (m9) ^_ 1> 



(4tt) b / 2 
2 

£ 



(19) 



/ 2 1 2 1 2\ 

2 (m x +m 2 +m 3 ) 



-(m 2 + m\ + mf) - Li(mi,m 2 , m 3 ) 



Li2(m\, rri2, ma) — 6£i(mi, m 2 , 7713) + (m 2 + m 3 — m 1 )ln(m 2 )ln(m 3 ) + (m 3 + mj — m2)lxi(m 3 )lxi(mi) 



+ (m 2 + ml - 77^)ln(mi)ln(m|) + £(mi,m 2 ,m 3 ) + (m\ + m\ + m 2 )[7 + C(2)] 



where we have defined 



ln(m 2 ) = In 



+ 7-ln(47r), 



Li(mi, TO2, 777,3) = m 2 ln(m 2 ) + m|ln(m|) + 777 2 ln(777 2 ), 

L2(77ll, 7772, 777 3 ) = 777 2 [ln(r77 2 )] 2 + 771?, [ln(?77 2 )] 2 + 777 2 [ln(?77 2 )] 2 , 

£ (mi,7772,m 3 ) = 4(2m 2 m 2 +2m 2 m 2 + 277737773 - m 4 - - m^) 1 / 2 L{0{) + L(9 2 ) + L(6 3 ) - ^ ln(2) 
In the latter expression, L(t) is the Lobachevsky function, defined as 



L(t) = f 
Jo 



dx In cos : 



and the angles are given by 



arctan 



(mj + m 2 + m§) - 2m 2 



{^lm\m\ + 2t?7 2 777 2 + 2m\m\ 



(20) 



(21) 

(22) 
(23) 

(24) 



(25) 



(26) 



These expressions, as well as the function £(7771,7772,7773), are valid for a positive argument of the square root. This 
depends on the value of the masses. For a negative value, one has to substitute 



£(mi, 777 2 , m 3 ) = 4(m 4 + m^ + m 4 - 2m 2 m^ - 2m 2 m 2 - 2m 2 .m 2 .) 1/2 [-M(-9i) + M(9 2 ) + M(9 3 )] , 



M(t) 



dx In sinhx, 



04 = coth 1 



(m\ + m\ + m 2 ) — 2m 2 



(m\ + m\ + m\ — 2m 2 m 2 — 2m 2 m 2 i — 2m 2 .777 2 ,) 1 / 2 



(27) 
(28) 

(29) 



In the following, we shall always keep the symbolic notation £(mi, 7772, 7773). This will avoid to check the sign of the 
argument of the square root. It is not necessary to do so in the pure </> 4 -case where we know that m^ > m-r, and, 
then, where the representation ( p4f ) is valid. In the present case, because we do not known the fixed points, it is not 
possible to specify the value of the photon mass vs. the transverse or longitudinal mass. We however mention that 
several simplifications arise when two masses are equal, or when one, or two masses, are vanishing. For two equal 
masses, say mi = 7772, the argument of the square root becomes m 3 (l — 4m 2 /777 2 ). This case arises when investigating 
the 0(iV)-theory. The ratio of the masses is the ratio Am\/m 2 L which is always smaller than unity. The relevant 
equations are then (24), (E5I) and The case of vanishing masses is a little bit problematic for the extraction of 



the power e° because of a cancellation of a diverging part of M(t) with a corresponding lm,- in (gfj). For this reason, 
we also give the following integrals: 
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(/i e ) 2 7(mi,0,0) 
(jM e fl(mi,mi,0) 



1 r(2-£)/2)r(3-D)r(D/2- l) 2 2 (m\ 



(47r) B T(£)/2) '"' ' 



2\ - £ 



1 r(2-D/2)r(l-D/2) 



(47T) 15 

2 



D-3 



2\ ~ £ 



'1 > 2 
/J 2 



(30) 
(31) 



( At £ ) 2 (4 7 r) 4 J(m 1 , m 2 , 0) = --(mf + to 2 ) 



-(to 2 + to 2 ) - Li(mi,m 2 ,0) 



7 2 (toi, m2, 0) — 6Li(mi, m-2, 0) + 2rn 1 ln(mj)ln(m 2 ) + [ln(mj — m 2 )] (m 1 — to 2 ) 



21n(TO 2 — TO 2 )ln(TO 2 )(TO 2 — to 2 ) + 2(to 2 — TO 2 )Li 2 



m; — m\ 



+ (m 2 + m 2 )[7 + C(2)] + — (m 2 -m 2 ) , (32) 



where Li 2 (z) = X)fcj_fV* 2 i s the dilogarithm. It is a simple exercice to check that, using mi — m 2 m (P^), we recover 
the e-expansion of (|3l|), while using 777.3 = 0, we recover the e-expansion of (|30|). In the previous determination of (|3^), 
we have assumed that toi > 777 2 . The case to 2 > m\ is obtained from the former by the replacement m\ <-> ni2. 
This only affects the term of order e°. Comparing Eq. J32] ) with Eq. (p0|), we can also extract the way £(mi, m2, 7773) 
diverges as one of the masses, say 7713, goes to zero. For m\ > m 2l we have 

lim (mi, TO2, m3) + (m 2 — TO2)[ m ( m i) — m ( m 2)] m ( m 3)l = 
m 3 — >o 



(m^ — m 2 ) < ln(m 1 )ln(m 2 ) + [lii(mi — m 2 )\ — 21n(mj — m 2 )ln(m 2 ) + 2Li 2 



ml 



(33) 



while the opposite case mi < m 2 is obtained by permuting toi and m 2 in this equation. From (33) we can also 
determine the case of two simultaneously vanishing masses: 

(34) 



lim {^(mi, TO3, m3) + 2m 2 ln(TO 2 )ln(m3) — m 2 [ln(r77,3)] 2 } = m\ < ln(TO 2 )] 2 + — 
m 3 — >o ' [_ 3 

With all these definitions, we are now armed to compute the two-loop diagrams. They are given by 



00 = iv { [{n - i)j{mT) + j(mi)]2 + 2 [(iv - i} j(wt)2 + j(wl)2] } 



(35) 



OQ = \v*{P - l)J(m,) [(N ~ l)J{m T ) + J(m L )} 



(36) 



o-(£) 



— ) $ [{N - l)I(m T ,m T ,m L ) + H(m L ,m L ,m L )] . 



(37) 



e 



= 4 



$ 2 < 7(m 7 , m 7 , m^) 



4m 7 



[ J(m 7 ) 2 (m|, - 2to 2 ) + 2m 2 J(m 7 ) J(m L ) + 2(m| - to 2 ) 2 7(to 7 , m L , 0) - m£7(m£, 0, 0)] I , (38) 



<~~> = eV^S (-^-2) [2J(to 7 )J(to t ) - J(m T ) 2 + 7(to 7 , m T , to t )(to 2 - 4to|)] 



+2 



(to 2 + m| - m|) 



J(?77 7 ) J(tol) + 



(to 2 - m| + ro|) 



J(to 7 ) J(tot) — J(tol)J(tot) 
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-I(m T ,m L ,0) 



{rrij. 



ml? 



+ I(m,j,mT,mL) 



m 2 L ) 2 + m 2 (m 2 — 2m\ 



2m 2 L ) 



(39) 



where, compared to Eq. (g) , we have included in the same diagram the different transverse and longitudinal compo- 
nents, and where the vertices and corresponding multiplicities have been once again specified. 



IV. EFFECTIVE POTENTIAL 



Collecting the results from the previous section, the renormalization coupling constants are tuned to cancel the 
poles in 1/e 2 and l/e. Working with the effective potential, as we do here, logarithms are appearing. Some of them 
have poles in e coefficients. These logarithms with pole coefficients have been named dangerous poles by Chung and 
Chung pgfl . The cancellation of these dangerous poles is a non trivial consequence of the renormalizability of the 
theory. We can see it as follow: Asking for the cancellation of one-loop poles, we identify straightforwardly 

.<> = ^ + ™' , (41) 

z™ = .</§• < 42 > 

where a factor l/(47r) 2 has been absorbed in the definition of h (see the expansion (FT 



(2) (2) 

The cancellation of the poles at the two-loop order gives renormalization coefficients Z m 2 , Z v which depend 

on ln(mi), ln(mT) and z9~ and a renormalization coefficient Z g 2 ^ which depends on ln(mi), ln(mT), ln(m 7 ) and 

, Z^' . We find that, asking for the independence with respect to these dangerous logarithms, fixes Z^P , . This 
is non trivial since there are more conditions than renormalization coefficients, and this is the mentioned consequence 
of the renormalization. Everything together, the cancellation of the poles at the two-loop order gives then 

Z( 3 = [9 2 (N + 5) - ISge 2 + 54e 4 ] - 1 [g 2 (N + 2) - 8ge 2 (N + 2) - 6e 4 (5iV + 1)] (43) 
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1 



g Zf> = — [g 3 (N + 8) 2 - 18g 2 e 2 {N + 8) + 108 5 e 4 (iV + 8) 2 + 108e 6 (iV + 18)] 



(j _ L 3 3 (5iV + 22) - 12 5 2 e 2 (iV + 5) - 18 5 e 4 (5A + 13) + 18e 6 (7iV + 90)] , (44) 

N 2 
Zi 2) =g-^[-l8e 2 + g(N + 2)]+ge 2 N-. (45) 

(46) 

where, as for the one-loop renormalization constants, a factor l/(47r) 2 has also been absorbed in the definition of h. 

With these one-loop and two-loop renormalization constants, all the poles disappear from the e-expansion of the 
theory Up to this order (two loops), and, in view of applications of this theory to phase transitions in three dimensions, 
the one-loop order term has to be developed to the order e. The two-loop effective potential in the Landau gauge can 
then be written as 



V = V (0) + _£_ [V (l,0) +£V (l,e) ] 



2 



V {2 \ (48) 



(47T) 

where we have specified that a factor (4-7r) 2 is absorbed in the definition of H, and with 

y(0) = i m 2$2 + 9 M e $ 4 + ^ (4Q) 

2 4! g^ £ 

0i,o) = ^!{(iV-l)m 4 .[-3 + 2E(m|)]+TO|(-3 + 2E(m|)]+m 4 [-5 + 6E(m2)]}, (50) 



V^ £) =-ILl((N- l)mi {21 - 181n(m 2 ) + 6[ln(m|)] 2 + tt 2 } + ml {21 - 181n(m 2 ) + 6[ln(m 2 )] 2 + tt 2 } 
96 V 

+ 3m 4 {9 - 10E(m|) + 6[hi(m 2 T )} 2 + tt 2 }) . (51) 

(52) 

The term V( 2 ) is much longer to write. For this reason, we give its expansion on the ln-bases. With the definition 



2 2 2 

V {2) =M- £ EEE y £UMm|)]iRm|)PRm 2 )] fe , (53) 

i=0 j=0 fc=0 



we have (see however the remark after these equations) 
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°'°> ~ 432 3 m4 



18m 2 [g 2 m 2 [(N - 1)(N + 5)m^ + lB>m 2 L m% - (2N + 13)m 4 L ] 
+ ge 2 {12(m\ - m|) 2 (m| + m\) - 6m 2 [(19 N - Z\)m% + 2Am 2 T m 2 L + 7mf ] 
- 18m 4 [(N - l)m 2 T + m|] + (197V - 36)m 7 } - 108e 4 m 2 (m| - m|) [(N - l)m 2 T + 3m|]^ 
+ 6 5 $ 2 {7 ff 2 m 4 [2(JV - l)m% + (N + 8)m|] + 18e 4 (7m| - 7m|m 2 + llm|m 4 + 54m 7 )} 
+ tt 2 (g 2 m* [2(JV - l)m| + (N + 8)m|] [3(m| - m 2 ) + ff $ 2 ] 
+ 18.ge 2 m 2 {(m| - m 2 T ) 2 (m 2 L + m|) + 3m 4 [(AT - l)m| + m|] - 9m 7 } 
+ 18e 4 {-9m 4 (mi - m 2 T ) [(N - l)m 2 T + 3m|] + g$ 2 (2m| - m£m 2 + 5m|m 4 + 10m 7 )} 

+ 216(47r) 4 5 e 2 [m 2 (m\ - m 2 T ) 2 (m L ,m T , 0) + e 2 $ 2 (m 2 - m|) 2 / (0) (m 7 , m L , 0) 

+ 54.ge 2 ^m 2 {(iV - 2)(m 2 - 4m|)m 2 £(m 7 , m T , m T ) 
+ 2 [mj. - 2m|(m 2 + m\) + (m 2 - m 2 L ) 2 ] £(m 7 , m T , m L )} 
+ e 2 ^ 2 (m i L — Am 2 L rn 2 + 12m 4 )£(m 7 , m 7 , mj,) 



+ 6 5 3 $ 2 m 4 [(AT - l)£(m T , m T , m L ) + 3£(m L , m L , m L )] 

(2) _ Tflrp / 2 2 2 ,/Ar, o\^2 , o^2l 



(54) 



Fl °>° = ~ 12^? I 9 {N ~ ^ m ^-< + ( N + + 2 9*1 



+ 6.ge 2 {3(m| - m|) 2 - m 2 [2m| + m|(7Af - 9)] + 5(AT - l)m 4 } - 54e 4 (N - l)m 2 (m| - m 2 T )j , (55) 

^■ 2 °>° = 72^ { g2{N ~ 1)m ' i~ m L( Qm T + + m 2 T [3m 2 (N + 3) + 4 3 <i> 2 ]} 

- 9 5 e 2 [-2 (mi - m|) 2 m| + m 2 (N - 2)(m 4 + 6m 4 .) - 6(2A^ - 3)m 4 m|] 

-162e 4 (A^-l)m 2 (m|-m|)m^, (56) 

= ~ 12^ (^Xt"^ + 5)m 2 + (2N + 7)m 2 + (TV + 8) 5 <f> 2 ] 
+ 6.ge 2 m 2 {3(m 2 - m 2 T ) 2 - m 2 [5m 2 L + 2m%] + 5m 4 } 

- 6e 4 {m 4 [27(mi - m\) - 16.g$ 2 ] + 3.g$ 2 m|(2m 2 - m|)}^ , (57) 



K% = J^\9 2 <\-^ N + 5 ) m l + 3( N + 8)^t + (N + 17).g$ 2 ] 
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+ 18 5 e 2 m 2 [(m 2 - m 2 T ) 2 + 3m 4 ] 

- 9e 4 {2m 4 [27(m| - m|) - 7 5 $ 2 ] + 5 $ 2 m|(4m 2 - 3m 2 )}^ , (58) 

2 

y ooi = 7Ar(-9e 2 ml$ 2 - 3m 2 [ml + m 4 - 2m 2 (m| + 4e 2 $ 2 )] 

+ 12m* [ml + (N - l)m 2 T - 6e 2 $ 2 ] - m^AN -9)}, (59) 

2 

F ° °' 2 = 8m^ [ _18m 7 + e2$2 (" m ' + 8m l m 7 - 22m 2 m 4 + 40m 7 )] , (60) 
VIZ = ^ {g(N - l)m|m 2 (3m| + 5 $ 2 ) 

[-(ml - m 2 T ) 2 (m 2 L + m 2 T ) + 3m 2 (m\ + m%) - 3m 4 (m| + m 2 .) + m^] |, (61) 



9e 



F iSi = — [-« - m r) 3 + 3m 2 (ml - m 4 ) - 3m 4 (m 2 - m 2 ) + (iV — l)m«] , (62) 



4mf, 



+ m 4 [m|(-3m| + 14e 2 $ 2 ) + 3m|] + 3m 7 (m| - m|) + m 7 j, (63) 

where the function I^(mi,rri2,0) denotes the e — > non-diverging piece of (|32|), without the coefficient (/U £ ) 2 - 
Some part of this function may be put in the other terms of the expansion (|53|), because they lead, for example, to 
ln(mi)ln(m2). We have not proceed in this way because, in the scalar QED case, it is not yet known if m 7 > tut- 
One has to study first the fixed points of the system. At the critical point, one knows that, in the 4 — e formalism with 
e-expansion, m — > 0. At the critical point, one would then need to compare e 2 to g/6. If m 7 > m^, one is allowed 
to use (32) with mi = m 7 , m 2 = m.^. In the opposite case m 7 < my, one has to use ( |32| ) with m 2 = m 7 , mi = mj. 
We also note that J"(mi, m2, 0) generate terms containing ln(m 2 — to?,). These terms can never be put on the basis 
defined in (p3|). 

To see the coherence of our effective potential, we shall look at the </> 4 -limit, which can be reached using e 2 — ► 0. 



limit 



In the limit e 2 — » 0, which gives the <fi 4 model, the renormalized coupling constants and the effective potential are 
much simplified. The renormalization constants can be rea dily read off Eqs.(|40|)-(^7|) with e 2 = 0. For the effective 
potential, it is necessary to clean the results given in ©-© and Taking the limit e 2 -> 0, we obtain 



V<® = im 2 $ 2 + 4 M e $ 4 + — , (64) 
2 4! gfi £ 

v (i,o) = ^|( iV -i) m 4 [-3 + 2bi(m 2 .)]+mi(-3 + 2hi(ml)]}, (65) 
8 

v (i,e) = _!Ll((N - l)mi {21 - 18m"(rr4) +6[ffi(m|)] 2 +tt 2 } + m£ {21 - 18ffi(m|) + 6[ffi(m|)] 2 + n 2 }) . (66) 
96 V / 

(67) 

The term is considerably simplified compared to the full QED model. It has the expansion 

2 2 

v(2) = ^"EE^H^JFNOF, (68) 

with the following coefficients: 

V$ = ^(6 {54mim| - 3(27V + 13)m£ + 7(N + 8)m 2 L g<f> 2 + (N - l)m 2 T [3 (AT + 5)m| + 14 5 i> 2 ] } 



10 



+ [2(N - l)m 2 T + (N + 8)ml] [3{<m% - ra\) + g<$> 2 ] tt 2 + 6g$ 2 [(N - l)£(m T , m T , m L ) + 3£(m L , m L , m L )]\ 



-^{m%(N - l)(N - 9) + mi(5N + 43) + m 2 L m 2 T (7N - 52) 



,2 m 2 



to t) [( n ~ l)£( m T,m T ,m L ) + 3£(m L ,m L ,m L )]} (69) 



V$ = "^(JV - %t [-mi + (JV + 3)m 2 T + 2g<S> 2 ] 

= -^(iV- l)m| [5m 2 + (iV-3)my] (70) 

^2,0 = ~^(JV - 1) {ml(Qm 2 T + ,g$ 2 ) - m|[3(iV + 3)m 2 + 4 5 0> 2 ]} 

= -^(^-l)K-3m|m 2 ,-(7V„l) m 4] (71) 

V$ = + 5)m 2 + (27V + 7)m 2 + (N + 8) 5 4> 2 ] 

= ^m£ [(N + l7)m 2 T - (2N + 19)mf ] (72) 
V$ = [-3(JV + 5)m 2 + 3(iV + 8)m 2 + (N + I7)g&] 

= |m|(4m! - 3m 2 ) (73) 

= ^(N-l)mi, (74) 

(75) 

where, in the second equality, we have replaced g& 2 by 3(m| — TOy). We have checked that Eqs. (]68|)-(|7^) reproduce 
the effective potential given in Q. The correct limit of the pure </> 4 -theory is then recovered. 



Remark 

In the previous subsection, we have obtained the effective potential of <j> -theory by taking the limit e 2 — ► for 
fixed m 2 , then taking the limit m 2 — > 0. This is the fastest way to obtain the correct limit, but it is quite cavalier to 
proceed so. We have checked, replacing everywhere m 2 by its value m 2 = e 2 $ 2 , that it was justified to do so: all the 
m 2 -terms in the denominator of Eqs. Q)-(^) are properly compensated to give the two-loop result of Eqs (|69|)-([74|). 



V. CONCLUSION 



In this paper, we have determined the effective potential of scalar QED in 4 — e dimensions in an e-expansion. For 
zero charge, we revover the well-known result for the pure 0(-/V)-symmetric 4 -theory. The full effective potential 
will be used in order to determine various amplitude ratios in the context of a scalar field in the presence of a gauge 
field. This comprises the superconducting case, for which N = 2. In this respect, it might be interesting to work in 
an arbitrary gauge and to show that these amplitude ratios are gauge independent. This can be done as in Ref. |p9| . 

We shall also, in a subsequent paper, recast our result in the scaling form of Widom which shows the critical 
behavior of the free energy explicitely. 

An other interesting work will be the application of variational perturbation theory to derive the various critical 
exponents in scalar QED, including the amplitude ratios along the lines of Ref. p7j . In 4 — e dimensions with 
e-expansion, this has not even been done for the simpler 0(7V)-symmetric </> 4 -model without gauge field. 
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